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Context. Proper motion measurements provide unique and powerful means to identify the driving sources of mass outflows, of 

particular importance in regions with complex star formation activity and deeply embedded protostars. They also provide the necessary 

kinematic information to study the dynamics of mass outflows, the interaction between outflows and the ambient medium, and the 

evolution of outflows with age of the driving sources. 

Aims. We aim to take a census of molecular hydrogen emission line objects (MHOs) in the p Ophiuchi molecular cloud and to make the 

first systematic proper motion measurements of these objects in this region. The driving sources are identified based on the measured 

proper motions, and the outflow properties are characterized. The relationship between outflow properties and the evolutionary stages 

of the driving sources are also investigated. 

Methods. Deep H? near-infrared imaging is performed to search for molecular hydrogen emission line objects. Multi-epoch data are 

used to derive the proper motions of the features of these objects, and the lengths and opening angles of the molecular hydrogen 

outflows. 

Results. Our imaging covers an area of ~ 0. 1 1 deg^ toward the L 1 688 core in the p Ophiuchi molecular cloud. In total, six new MHOs 

are discovered, 32 previously known MHOs are detected, and the proper motions for 86 features of the MHOs are measured. The 

proper motions lie in the range of 14 to 247 mas/yr, corresponding to transversal velocities of 8 to 140km/s with a median velocity of 

about 35 km/s. Based on morphology and proper motion measurements, 27 MHOs are ascribed to 21 driving sources. The molecular 

hydrogen outflows have a median length of ~0.04 pc and random orientations. We find no obvious correlation between H2 jet length, 

jet opening angle, and the evolutionary stage of the driving sources as defined by their spectral indices. We find that the fraction 

of protostars (23%) that drive molecular hydrogen outflows is similar to for Class II sources (15%). For most molecular hydrogen 

outflows, no obvious velocity variation along the outflow has been found. 

Conclusions. In Ophiuchus the frequency of occurrence of molecular hydrogen outflows has no strong dependency on the evolutionary 

stage of the driving source during the evolution from the protostellar stage to the classical T Tauri stage. 

Key words, stars: formation - stars: winds, outflows - ISM: jets and outflows - infrared: ISM - shock waves 
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1. Introduction 

Mass outflow plays an important role during various stages of 



star formation (ShuetaL 1987 Arce et al. 2007 Bally et al. 



^ |2007| l. It is believed that mass outflow is a vital channel for the 



transfer of excess angular momentum from the molecular cores 
to the ambient interstellar medium (Shang et al.'2007). Outflows 
are also good tracers of the young stellar objects (YSOs), in 
particular those deeply embedded in the molecular cloud cores. 
They can be studied in various regions of the electromagnetic 
spectrum. In the optical, Herbig-Haro (HH) objects are the man- 
ifestation of shock-ionized mass outflows, and trace either ma- 
terial ejected from the protostars directly or the shocked in- 
terstellar medium ( Reipurt h & Bally|l2001t |Wang et alT][2004 
Walawe nder et al.||2005[ |Wang et al.||2005[]Wang & Henning 



2006, 2009 1. The typical velocity of HH objects is up to several 



* Based on observations made with ESO New Technology Telescope 
at La Silla under programme ID 079. C- 0717(B), and on data obtained 
from the ESO Science Archive Facility. 



hundred km s '. At millimeter wavelengths, more than 400 CO 
molecular outflows have been detected ( |Wu et al.|2004| l. These 
CO outflows probe the swept-up or entrained medium in the out- 
skirts away from t he central jets, a nd have typical velocities of a 
few to ten km s"' ( |Bachiller|l996| . 

In the near-infrared, molecular hydrogen emission lines, in 
particular the v = 1-0 S(l) transition at the wavelength of 
2.12 yum, are powerful tracers of shock-excitation. Davis et al. 
( |2010) strictly defined molecular hydrogen emission-line ob- 
jects (MHOs) as the near-infrared manifestations of outflows, 
and presented a comprehensive catalog of around 1000 MHOs. 
Deep, wide-field near-infrared imaging with the combination 
of narrow-band H2 filter and the corresponding filter for con- 
tinuum emission such as a Ks filter has become a very effi- 
cient method to detect and identify MHOs (jStanke et al.]|2002 



•Rhanzad yan et al.|2004[ [Davis et al.|20"08||Froebrich et al.|201 l| l 
Morphology, proper motion, as well as association with an HH 
object or CO outflow can all provide information concerning the 
location of the driving source. Once the relationship between the 
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MHOs and their driving sources has been established, we can 
study the distribution of YSOs, the interaction of YSOs with the 
ambient interstellar medium, the star formation efficiency, and 
the evolution of the entire star forming region by analyzing the 
statistical characteristics of the MHOs. MHOs in several nearby 
star forming regions have been investigated using deep, sub- 
arcsecond-resolution, near-infrared surveys ( [Stanke et al.|2002| 
[Davis etal..2008, 2009) . 

The p Ophiuchi molecular cloud is the nearest star form- 
ing region at a distance of about 120 pc ( Lombardi et al.|20W 



Loinard et al.|2008]l. Lynds ( 1962 1 identified several dark nebu 



lae in Ophiuchus by studying red and blue prints of the National 
Geographic -Palomar Observatory Sky Atlas. The large-scale 
structure of the p Ophiuchi molecular cloud h as been revealed 



by extensive molecular '^CO line mapping dLoren 



1989a 



bl 



Surveys for young stars in Ophiuchus include] Casanova et al.| 



(fT995l l [Kamata et al.| ^\WI), [Grosso et al. (2000), Imanishi 
let al.| ( |2002li, |Ozawa et al.| (|2005| l, a nd |Pillitteri et al.,(,2010j 
in X-rays; [Wilking et al.|(|1987a|l and IWilking et al.| (|2005[ ) in 



the visual; [Greene & Young' ('1992a' 



|2004[ [Gomez et al.|2003||Khanzadyan et al.|2004| i, which have 
large uncertainties. The kinematic information of the MHOs, in- 
cluding proper motions and radial velocities, will offer reliable 
evidence for the identification of the driving sources in this re- 
gion. 

This work is part of a project aimed at understanding the dy- 
namical characteristics of outflows in Ophiuchus. The objectives 
of the project are to measure the radial velocities of jets with 
VLT/CRIRES, and to derive the proper motions of jets based on 
multi-epoch astrometry with NTT/Sofl. This then in turn should 
enable us to derive the flow orientation, trace the driving sources 
of the jets, and based on the 3D velocity information and spec- 
tral line features to compare observed jet properties with models 
of shock physics. As a first step, we here present the proper mo- 
tion measurements of MHO features in the p Ophiuchi molecular 
cloud. The radial velocity measurements of MHO features with 
VLT/CRIRES in Ophiuchus will be published in another paper 
(Zhang et al. 2013, in preparation). 



Alves de Oliveira & Casali (2008) in the near-infrared; [Young[ 



Allen etai.[ ([2 002^, and 2. Observations and data reduction 



et al.[( |1986a), Wilking etaL[([20()T i, [Padgett et a l. (2008), Evans 
eraL|p009) and|McClure et al. (|2010|) in the mid- and far 



[et al 



infrared; and[Motte et al. ( 1998 1, Johnstone et al. 
^rar('2006^,IAndrews & Winiams'f2007), and 



(|2000li,|Stanke 



Andrews et al. 



2010) at sub-millimeter and millimeter wavelengths. By com- 
bining the results of the previous studies, |Wilking et al.[ ( [2008] l 
compiled a list of 316 young stars in L1688, the densest core of 
Ophiuchus. Using observations by the Spitzer space telescope, 
[Evans et ar] ( [2009 1 derived the most complete and unbiased YSO 



sample for the overall p Ophiuchi star forming region. Evans 
et al. ( 2009|l investigated the overall evolution and estimated for 



the p Ophiuchi cloud a star formation rate (SFR) per unit area 
of about 2.3 Mq Myr ' pc"^. This is lower than the SFR in the 
Serpens cloud (3.2 M© Myr ' pc"^) but h igher than the SF R in 
the Perseus cloud (1.3 M© Myr"' pc"^) ( ,Evans et al. 2009 1. 

Outflows in Ophiuchus have been investigated in some stud- 
ies at different wavelengths. Close to 50 Herbig-Haro (HH) 
objects ( Wilking et al. 19971 [ Gomez et aL][1998[ [Wu et al. 
[2002; Phelps & Barsony 2004) and more than 15 high- velocity 



CO molecular outflows (|Andre et al. 



Bontemps et al. 



Nakamura et al. 



1996) [Wu et al.[|2004[ [Bussmann et al. 



[T993i [Dent et al. 



1995^ 



2007^ 



2011 1 have been identified. Using mid-infrared 



Spitzer/IRAC observations, [Zhang & Wang[ ( |2009] l detected 44 



extended green object (EGO) outflow features. Based on near 
infrared data, Davis et al. (2010|) list ~47 MHOs consisting of 
~119 H2 2.12 jum emission features in the p Ophiuchi molecu- 
lar cloud. While observations with small angular coverage de- 
tected some individual MHOs ( [Grosso et al. 2001^ Ybarra et aT] 
[2006| l, the majority of MHOs in Ophiuchus have been detected 
by two unbiased surveys ("Gomez etal.|2003t|Khanzadyan et al 



2004^. We observed a subset of these MHOs, aiming at the sec- 



ond epoch images to measure the proper motions (PMs) of the 
MHO features. 

Compared with outflow morphology, motions of outflow fea- 
tures are more robust evidence for the identification of the driv- 
ing sources of outflows. The directions of motion vectors indi- 
cate the location of the driving sources of the outflow features. 
L1688, the densest region in Ophiuchus, includes hundreds of 
young stars and several tens of outflow features within an area 
of less than one square degree. The relationship between the 
young stars and the outflow features is very complex in this re- 
gion. Previous studies associated the outflow features with the 
young stars just based on the morphologies ( |Phelps & Barsony| 



2. 1 . Near-IR imaging 

The observations were conducted in service mode between July 
and September 2007 (079.C- 0717(B), PI: M.D. Smith) with the 
ESO NTT near-infrared spectrograph and imaging camera SofI 
( Moorwood et al.|l998 1, covering in total an area of ~ 0. 1 1 deg^ . 
Figure [11 shows the coverage of the observations with the gray 
filled boxes in the left panel. SofI is equipped with a Hawaii 
HgCdTe 1024x1024 array, yielding a field of view of -4:9x4^9 
with a plate scale of 0'.'288 pixel ' in its large field (LF) imag- 
ing mode. In total, 16 fields were observed in K, and H2 band. 
The broadband K, filter has a central wavelength of 2.162 /urn 
and a bandwidth of 0.275 fim. The narrowband H2 filter has 
a central wavelength of 2.124 /urn and a bandwidth of 0.028 
//m. Using dithered observations, each field was typically im- 
aged seven times with the individual exposure time of 18s in Kj 
band and of 180s or 240s in H2 band. Table[T]lists the center co- 
ordinates, the total integration time, and the average airmass for 
each field. 



Table 1. Observing log 



Center of field 


Exposui 


e time 


Airmass 


RA 


Dec. 


H2 


Ks 


H2 


Ks 


(J2000) 


(J2000) 


(s) 


(s) 






16 26 16.0 


-24 24 30 


1680 


126 


1.02 


1.01 


16 26 21.0 


-24 27 55 


1680 


126 


1.05 


1.04 


16 26 33.0 


-24 25 20 


1680 


126 


1.10 


1.09 


16 26 53.9 


-24 37 33 


1260 


126 


1.18 


1.17 


16 27 03.3 


-24 47 03 


1260 


126 


1.27 


1.26 


16 27 08.8 


-24 33 43 


1680 


144 


1.49 


1.47 


16 27 11.8 


-24 43 59 


1260 


126 


1.75 


1.72 


16 27 22.4 


-24 49 13 


1260 


126 


1.08 


1.07 


16 27 24.3 


-24 28 17 


1680 


198 


1.14 


1.12 


16 27 34.1 


-24 39 18 


1260 


126 


1.64 


1.61 


16 27 43.0 


-24 40 47 


1260 


126 


1.24 


1.23 


16 27 45.5 


-24 45 42 


1260 


126 


1.37 


1.35 


16 27 47.7 


-24 3150 


1260 


126 


1.89 


1.85 


16 28 01.6 


-24 3102 


1260 


126 


1.53 


1.51 


16 28 10.2 


-24 16 01 


1260 


126 


1.92 


1.88 


16 28 20.3 


-24 37 07 


1260 


126 


1.37 


1.36 



We processed the data using the IRAF packages. After sub- 
tracting a median dark frame, each image was divided by a nor- 
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malized flat-field for each filter. We then corrected the illumi- 
nation effect and masked the bad pixels using the calibration 
data available at ESO's websita^ and then subtracted the sky 
frames, which had been generated from a median-filtered set of 
the flattened frames in Ks and H2 bands. For the SofI imaging 
mode, a bright source on the array produces a "ghost" that affects 
all the lines where the source is and all the corresponding lines 
in the other half of the detector, a instrumental feature called 
"Inter-quadrant Row Cross Talk" (details can be found in SofI 
User's ManuaQ. We also corrected this instrumental effect with 
the IRAF script from SofI websitaj In order to get the accurate 
astrometry, we also removed the geometrical distortion on the 
dithered frames using the geometric distortion solution. Finally 
we shifted and combined the dithered frames in each filter for 
each field after excluding several images with bad quality. 



2.2. Proper motion measurements 

Archival SofI images ((PID: 65.1- 0576(A) and 67.C- 0284(A)) 
obtained in July 2000 and August 2001, and published by 



Gomez et al. ( j2003) ; |Khanzadyan et al. ( 2004 1 served as the first 
epoch data. The newly-acquired SofI images provided the sec- 
ond epoch to measure the proper motions of the H2 emission- 
line features in MHOs. We reprocessed the data using an IDL 
pipeline which was improved by |Gennaro et al. (2012 1 accord 



2.1 



ing to the process of data reduction described in Section 
Figure [T] shows the field coverage. The gray filled boxes show 
the coverage of the Soft (2007, PID: 079.C- 0717(B)) observa- 
tions, while the red dashed and green solid polygons outline the 
coverage of the Soft 2000 (PID: 65.1- 0576(A)) and 2001 (PID: 
67. C- 0284(A)) observations, respectively. 

In order to obtain accurate proper motion measurements, we 
firstly corrected for seeing. For each image in the image pair 
(first epoch image and corresponding second epoch image), the 
average full width at half maximum (FWHM) of all detected 
stars is used to represent the seeing condition. Then the im- 
age, which has the smallest average FWHM is convolved with 
a Gaussian core so that two images in each image pair have the 
same FWHM. After convolution, we align and resample the im- 
ages in order to match the image pairs pixel to pixel. The as- 
trometric calibration was performed using the IRAF packages. 
Using stars in common between the different data sets, we reg- 
istered the first epoch SofI images to the 2MASS Point Source 
Catalog, and then registered each second epoch SofI image to the 
corresponding first epoch image. The root-mean-square (rms) 
values for the registration residuals of all fields were below 
0.1" (-1/3 pixel). 

The proper motions of the MHO features are measured with 
a method similar to the one suggested by McGroarty et al. ( 2007| l 
and |Davis et al.| ( |2009[ ). For each MHO feature, we defined man- 
ually a polygon aperture based on the surface brightness distribu- 
tion to estimate the flux and morphology of the feature. Then the 
first epoch image was shifted relative to the second epoch image 
in steps of 0.25 pixels about an 8 x 8 pixel grid. At each position 
the flux of the features in the polygon are used to calculate the 
cross-correlation coefficient. This process yields a set of correla- 
tion coefficients as a function of shifts (ccf - f{shiftx, shifty)). 

' http://www.eso.org/sci/facilities/lasilla/instruments/sofi/tools 
/reduction/index.html 

^ http://www.eso.org/sci/facilities/lasilla/instruments/sofi/doc/manual 
/sofiman_2pl0.pdf 

^ http://www.eso.org/lasilla/instruments/sofi/tools/reduction 
/sofi_scripts/ 



Obviously, the shift corresponding to the maximum correlation 
coefficient can be obtained as the shift of each feature between 
two epochs. Note that we determined the position of the maxi- 
mum correlation coefficient by calculating an intensity weighted 
center using pixels in the correlation coefficient image around 
the peak value. Finally, we subtracted the apparent proper mo- 
tions of the adjacent background field stars from this shift in 
order to minimize any systematic effects such as the problems 
with alignment, resampling, or astrometry. As for the method 
described in detail by Davis et al. ( 2009 1, the errors of PMs of 
flow features are mainly from the random motions of the field 
stars. 

The discussion above assumes that the brightnesses and mor- 
phologies of the MHO features do not change between the first 
and second-epoch observations: for many features this is not the 
case. In fact, many factors such as signal-noise-ratio, shape, flux 
variability of the MHO features, and the gradient in the back- 
ground emission also contribute to the errors of the PM mea- 
surements, but these factors are very hard to quantify. Therefore, 
we checked every feature and introduced a parameter called 
"Status" to estimate the reliability of the PM calculations. A PM 
with status is reliable, while a PM with status 1 or 2 is uncer- 
tain. 

The proper motions of MHO features obtained through the 
above process are relative to the background field stars. 'MakarovJ 
(2007) investigated 58 probable members of the association 
of pre-main-sequence stars around the filamentary p Ophiuchi 
cloud using astrometric proper motions from the UCAC-2 cat- 
alog and the convergent point method. They found that the se- 
lected probable members have similar proper motions with an 
average of roughly (jj.aCos6,ius) = (-11.4,-24.5) mas yr"'. The 
bulk of the standard proper motio n errors are in the range 1-12 



mas yr per coordinate direction. Mamajek (2008 1 also inves- 



tigated 53 candidate members in L1688 based on the Tycho-2, 
UCAC2 and SPM2.0 catalogs. They calculated the mean proper 
motion for members in the group with several estimators which 
are fairly insensitive to outlying points and finally obtained the 
mean proper motion of {jJaCos6,fj.s) - (-10, -27) mas yr"' with 
the total uncertainty of (2, 2) mas yr"'. We adopted the mean 
proper motion from Mamajek (2008) to represent the system- 
atic proper motion of the members in L1688. We subtracted this 
systematic proper motion from all the PM measurements of the 
MHO features with the aim of getting the proper motions of the 
MHO features with respect to the p Ophiuchi molecular cloud. A 
complete list of proper motions of MHO features in p Ophiuchi 
is provided in the Appendix. 

2.3. Flux estimation of MHO features 

Areal photometry was performed to estimate the flux of MHO 
features. We firstly defined a polygon aperture based on the mor- 
phology and surface brightness distribution of each MHO fea- 
ture. A 10 pixel wide region of pixels surrounding the polygon 
aperture is used to estimate the sky background. The polygon 
apertures and the surrounding annulus are applied to the MHO 
features on the continuum-subtracted H2 images to estimate the 
fluxes of H2 1-0 S(l) emission line. 

The flux calibration was obtained via the 2MASS catalog. 
We used Sextractor ( |Bertin & Arnouts 19961 to do the pho- 
tometry for the point sources on the H2 images. The zero-point 
flux was derived through the comparison of the integrated counts 
with the cataloged flux of the point sources. We applied the zero- 
point flux to the integrated counts of the MHO features with the 
aim to convert the counts of MHO features to the flux unit. 
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Fig. 1. Coverage of our data. In the left panel, the underlying image is 2MASS image in Ks-band. The gray filled boxes show 
the coverage of Sofl (2007, PID: 079. C- 0717(B)) observations. The red dashed and green solid polygons outline the coverage 
of the Sofl 2000 (PID: 65.1- 576(A)) and 2001 (PID: 67.C- 0284(A)) observations, respectively. The YSOs identified by Spitzer 
photometry (Evans et al. 2009 1 are marked with filled circles. In the right panel, the gray contours represent the continuum emission 
at 11 20 yum from the COMPLETE project ( jRidge et al.|2006p with the contour levels from 80 to 2000 mJ y/beam and a step of 100 
mJy/beam. The green filled triangles show the positions of the 1 100 jum millimeter sources identified by |Young et al.| ( i2006) . The 
newly detected MHO features are marked with black open squares and the previously known MHO features are labeled with red 
pluses. The YSOs identified by Spitzer photometry are marked with blue circles. 



Figure. l2] shows the distribution of the fluxes of the MHO 
features. The median of the fluxes of MHO features is about 
6x10 '** W mT^. The minimum of the fluxes we detected is about 
6xlO-'''Wm-2. 
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Fig. 2. Histogram of integrated 1-0 S(l) fluxes for the detected 
MHO features. 



3. Results 

3.1. The detected MHOs in Optiiuctius 

We performed our survey towards the L1688 core, which is the 
densest region in the p Ophiuchi molecular cloud, covering an 
area of about 0.1 deg^. We detected six new MHOs and 32 
known MHOs, with in total 107 MHO features. Table|2]lists the 
position and morphology of six newly detected MHOs and figure 
[T]shows the distribution of all 107 MHO features. 

Note the diff'erence between "MHO" and "MHO feature". 
An MHO feature is an H2 emission feature which can not be 
resolved to more sub-structures in H2 images. An MHO consists 
of one or more MHO features. Therefore, we detected 38 MHOs 
which consist of 107 MHO features and the PM measurements 
and photometry mentioned above are calculated for the MHO 
features. 

MHO 2149, MHO 2154, and MHO 2155 are faint feahires 
which are only detected in our deep Sofl (2007) images. MHO 
2151-2153 are visible in two epoch images. We reprocessed the 
Sofl (2001) data and detected MH021 51-2153, which were not 
identified by' Khanzadyan et aL] ( |2004] ) in the first epoch images. 
Therefore, we also obtained the proper motions of MHO 2151- 
2153. The details of all detected MHOs in Ophiuchus (including 
names, positions, proper motions and the notes on individual ob- 
jects) can be found in the Appendix. 



3.2. The proper motions of MHO features 

We obtained the proper motion measurements for 86 MHO fea- 
tures in Ophiuchus. Table A.l in the Appendix lists the proper 
motion of each MHO feature. Here we present some statistics of 
these PM measurements. 
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Table 2. Newly detected MHOs in Ophiuchus 



MHO names" 


RA 


Dec. 


Comments 




(J2000) 


(J2000) 




MHO 2149 


+ 16 26 17.9 


-24 27 29.4 


faint knot 


MHO 2151 


+ 16 26 21.2 


-24 25 00.3 


three faint knots in a diffuse nebula 


MHO 2152 


+ 16 27 13.3 


-24 26 31.9 


knot with a faint diffuse nebula 


MHO 2153 


+ 16 27 20.1 


-24 26 25.5 


knot 


MHO 2154 


+ 16 28 16.1 


-24 34 58.2 


faint diffuse nebula 


MHO 2155 


+ 16 27 24.6 


-24 29 12.0 


faint elongated knot 



Notes. '"' These names are provided by C.J. Davis and will be included in the online MHO catalog ihttp://www.astro.ljmu.ac.uk/MHCat/J. Please 
note that there is no MHO 2150 in this table. Actually, MHO 2150 has been detected by Caratti o Garatti et al. (2006', see Fig. 12 in their paper), 
but somehow not been included in the MHO catalog before. We thus associated this object with MHO 2150 and included it in Table. [a!T] as a 
previously known MHO. 
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Fig. 3. Top: distribution of tangential velocities of MHO features 
(in 10 km s"' bins); middle: the distribution of fiacosd (in 30 
mas yr"' bins); bottom: the distribution of /i^ (in 30 mas yr"' 
bins). The open columns are for all PM measurements of MHO 
features in Ophiuchus and the gray filled columns are for the 
PM measurements with status of '0' (reliable measurements) of 
MHO features. The red dashed line in top panel shows the fit- 
ting with the exponential function for the MHO features with 
the reliable tangential velocities greater than 40 km s"'. 



We found that the proper motions of MHO features are in 
the range of 14 to 247 mas yr"', corresponding to transver- 
sal velocities of 8 to 140 km s ' if assuming that all of MHO 



features are located at a distance of 119 pc (Lombardi et al. 
12008 ). Figure. [3] shows the distribution of tangential velocities 
(top panel), jJaCosS (middle panel) and fig (bottom panel) for all 
PM measurements of MHO features with the open columns and 
for reliable PM measurements (Status "0") with the gray filled 
columns. The median of the distribution of tangential veloci- 
ties is 34.5 km s"' for all PM measurements of MHO features 
and 34.7 km s"' for the reliable PM measurements. The peak of 
the distribution of tangential velocities is in the range of 30-40 
km s ' with the variation in the bin size. Moreover, the MHO 
features with tangential velocities > 40 km s"' show an expo- 
nentially decaying frequency. Actually, we fitted the tangential 
velocity distribution of reliable PM measurements (marked with 
status 0) with the value of > 40 km s"' using the power-law func- 
tion and the exponential function, individually. The fitting with 
the exponential function returned the smaller chi-square value. 
Therefore, we think that the distribution of tangential velocities 
with the value > 40 km s"' is closer to the exponential distribu- 
tion. The top panel of Figure. l3] shows the fitting result with the 
red dashed curve. The counts of MHO features is related to the 
tangential velocities of MHO features in the way of Counts oc 

g(-0.071±0.011)xiWoc!7j(*:m/.s) 

Davis et al. ( 2009|l measured proper motions for 147 H2 fea 



tures in Orion A. The median of the distribution of tangential 
velocities is 70.8 km s"' for all PM measurements of MHO fea- 
tures and 75.3 km s ' for the reliable PM measurements of MHO 
features in Orion. The peak of the distribution of reliable PM 
measurements of these MHO features in Orion is in the range of 
60-80 km s"'. The distribution of tangential velocities with the 
value > 70 km s ' also shows an exponential distribution in the 

way of Counts oc Qi-0-028±0.004)xvelocity(km/ s) ^ 

We also note that there are drops in number at low veloc- 
ity in the distributions of tangential velocities of MHO features 
in Ophiuchus and in Orion. These drops should result from 
the incompleteness of the PM measurements. For the low sur- 
face brightness MHO features and the slow moving MHO fea- 
tures, we cannot obtain reliable proper motion measurements. 
Therefore, these drops are due to a bias. More sensitive observa- 
tions with higher angular resolution are required to construct the 
correct distribution in the low velocity range. 

3.3. The driving sources of MHOs in Ophiuchus 

The most probable driving sources of MHOs are identified based 
on the outflow morphology and proper motion analysis. In to- 
tal, we associated 27 MHOs with 21 YSOs identified by 



etaL| ( [2009l ). Table|3]lists these 27 MHOs and their likely driving 



Evans 



sources. The driving sources of the remaining eleven MHOs can 
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not be identified because of the uncertainties in their proper mo- 
tions. Figure HI shows the distribution of the 21 driving sources 
of MHOs in Ophiuchus. 

Note the difference between the number of MHOs (27) and 
the number of molecular hydrogen outflows (21). A molecular 
hydrogen outflow may consist of two or more MHOs that are 
associated with the same driving source. 

Previous studies (Gomez et al. 2003 Khanzadyan et al. 



2004 1 also associated MHOs in Ophiuchus with some young 



stars, but their identifications are only based on MHO morphol- 
ogy and there are large uncertainties. Our proper motion mea- 
surements supply more robust evidence for the driving source 
identifications of MHOs in Ophiuchus, especially for the com- 
plicated regions that have a large number of young stars around 
the MHOs. For example, in the region with VLA1623 (see 
Figure. |A.3| l, there are several tens of MHO features and tens 
of young stars. Previous studies (Dent et al. 1995; Gomez et al. 



2003) suggested that most of the MHO features in this region 



should belong to the outflow which is driven by VLA 1623 based 
on the association between the MHO features and the CO out- 
flow driven by VLA 1623. But they cannot point out which fea- 
tures are exactly excited by VLA 1623 and which are not, espe- 
cially for the features near the young star GSS 30. We derived 
reliable PM measurements for most of MHO features in the re- 
gion of VLA 1623 (also see Figure. |A.3[ ). We found that most 
of the PM vectors trace back to VLA 1623, which confirms that 
these MHO features are driven by VLA 1623. We also found that 
there is no obvious evidence for the assumption that some MHO 
features in this region are excited by GSS 30. 

Among 24 previously known MHOs which are associated 
with the driving sources in our paper based on the morphologies 
and the PM measurements, there are 13 MHOs (~54%) whose 
driving sources are in agreement with the previous identifica- 
tions while only 5 MHOs (~21%) are associated with the differ- 
ent driving sources from the previous identifications. Note that 
there are 4 MHOs (~17%) which were associated with several 
possible driving source candidates in previous studies. Now we 
can provide a robust identification for their driving sources based 
on our PM measurements. Note that there are two previously 
known MHOs (MHO 2103 and MHO 2104) for which we did 
not obtain the PM measurements because of the lack of the first 
epoch SofI data. 

|Nakamur a et al.| p011| l have published the results of CO 
mapping observations toward the L1688 core in Ophiuchus. 
They identified five CO molecular outflows based on the CO (J 
- 3-2) and CO (J = 1-0) emission maps. We also checked for the 
associations between MHOs and CO outflows. Table [3] lists the 
results. We found that two of the driving sources of MHOs also 
drive the CO outflows. The MHOs driven by VLA 1623-243 ai-e 
well associated with the CO outflows detected in the CO (J - 
3-2) emission map. Many MHO features that are driven by VLA 
1623-243 correspond to the red-shifted lobe of the CO outflow. 
[Nakamura et al.|p01 l| l also detected a weak CO outflow on their 
CO (J = 1-0) map and they associated this CO outflow with the 
YSOIRS44 . This CO outflow was firstly identified by |Terebey| 



et al. ( 1989 1 by means of CO (J = 1-0) interferometric maps us- 



ing the Owens Valley Millimeter Interferometer We have com- 
pared the position and the morphology of their CO outflow (see 
Terebey et al.|1989 , Fig. 2b) with our MHOs and we suggest to 
associate this CO outflow with YLW 16 - which also drives two 
MHOs - rather than IRS 44. We also detected an MHO, MHO 
2153, corresponding to the blue-shifted lobe of the CO outflow 
that is driven by Elias 2-32. However, according to the proper 
motions we obtained for MHO 2153, we associate MHO 2153 



with [GY92] 239 rather flian with Elias 2-32. For the CO outflow 
driven by Elias 2-29, we can not detect any MHO counterparts as 
our H2 SofI (2007) map does not cover the region around Elias 
2-29. 

Johnstone et al. ( 2000| l have published the results from a 
survey of the central 700 arcmin- region of the p Ophiuchi 
molecular cloud at 850 /im using the submillimeter common- 
user bolometer array (SCUBA) on the James Clerk Maxwell 
Telescope and they identified 55 cores on their 850 /zm map. 

((2006^ 



Moreover, Young et al. 



also published the results 
from large-scale millimeter continuum map at 1.1 mm of the 
p Ophiuchi molecular cloud with Bolocam on the Caltech 
Submillimeter Observatory (CSO) and they detected 44 definite 
sources. [Evans et aL] ( |2009[ l correlated these dust cores with their 
YSO sample in order to get the complete spectral energy distri- 
butions (SEDs). Among our 21 flow sources, four sources (19%) 
have 850 /vm core counterparts and four sources (19%) have 1.1 
mm core counterparts. In fact, there are about 184 YSOs in the 
L1688 core region, of which 23 YSO (~12.5%) are associated 
with 850 jum dust cores, and 22 YSOs (~ 12%) have 1.1 mm 
core counterparts. YSOs driving molecular hydrogen outflows 
tend to have a higher degree of association with millimeter/sub- 
millimeter cores than YSOs not driving molecular hydrogen out- 
flows. 



3.4. Outflow statistics 

Based on the positions and proper motions of the MHO features, 
we calculated the jet length (L), jet opening angle (6), jet position 
angle (PA), variation in the velocity along the outflow (VV), and 
dynamic ages (DA) for the molecular hydrogen outflows. We 
defined some parameters (L, 0, PA) according to the suggestions 
from |Davis eTaL] ( |2009l l. 

Jet opening angle (0) is measured from a cone with the small- 
est vertex angle centered on the driving source that includes all 
MHO features in the outflow. The opening angle can not be cal- 
culated if there is only one feature in the outflow. Jet position 
angle (PA) is measured east to north as computed from the bisec- 
tor of the opening angle. We selected the angle smaller than 180 
degrees as the position angle for the bipolar outflows. Jet length 
(L) is defined as the distance from the MHO feature to the source 
projected to the bisector of the opening angle in cases where 
more than one feature is identified in the outflow. In case where 
there is only one feature identified in the outflow, this length is 
just the distance from the source to the feature. Variation in ve- 
locity (VV) along the outflow is calculated through fitting the 
tangential velocity of MHO features vs. the distance from MHO 
features to the driving source using a linear least square method. 
Note that we calculate the VV of each lobe for the bipolar out- 
flows. VV equals 1 means that the MHO features will speed up 
by 1 km s"' after they propagate 100 AU along the outflow if 
assuming that the outflows are located at the distance of 119 
pc. The dynamic age (DA) of the outflow is the maximum dy- 
namic age of the features with proper motions. Table l3] lists 
these parameters for each outflow. Note that these parameters 
are only calculated with MHOs (not considering the information 
provided by HH objects and CO outflows) and uncorrected for 
inclination to the line of sight. 

As our data cover only part of Ophiuchus, we selected the 
region outlined in Figure|4]with a black solid polygon according 
to the coverage of unbiased surveys conducted by Khanzadyan] 
et al. ( 20^, [Gomez et al.| ( |2003 1 and the coverage of our Sofl 
(2007) data. In this region, our MHO sample is complete. This 
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area covers the majority of the L1688 core, and our follow-up 
statistics are all restricted to this area. 

In this area, there are 123 YSOs ( |Evans et aL][2009 ll, in- 
cluding 52 protostars (-0.5<q'<3.0) and 89 disk-excess sources 
(-2.5<Q'<0.0). Of these 52 protostars, 10 (19%) sources drive 
outflows. This fraction is lower than the fraction (~30%) of pro- 
tostars with outflows in Orion A ( |Davis et al.|2009| l. [Davis et al.| 
( 2009| l identified 116 jets in Orion A and most of their driving 
sources are protostars. However, of our 2 1 sources with outflows, 
10 are protostars and 12 are disc-excess sources. The criterion of 
the classification for protostars and disk-excess sources is from 



protostars are born, with the associated outflows aligned par- 
allel with the local field, and perpendicular to the chains of 



cores ([Mouschovias 1976 Banerjee & Pudritz|2006 Davis et al. 



Davis et al. (2009) and note that there is an overlap between the 



range of a for protostars and for disc-excess sources. 

Ophiuchus is located considerably closer to us (119pc, 
Lombardi et al.|2008| l than Perseus (-300 pc, [Bally et al.|2008] l 



and Orion A (-420 pc, [Menten et al.|2007[ ). Of our 21 molecular 
hydrogen outflows, 12 (57%) exceed one arcminute in length (at 
a distance of 119pc, 1' is equivalent to 0.035 pc) and 4 (19%) 
exceed three arcminutes (corresponding to 0. 1 pc at a distance 
of 119 pc). A histogram (top panel in figure B} shows the distri- 
bution of jet lengths. Shorter jets are more common than longer 
jets. The median of jet lengths is 1.1', corresponding to the phys- 
ical scale of about 0.04 pc in Ophiuchus. 

Davis et aL] p009) identified 116 jets in Orion A and the 



lengths of these 116 jets show an exponential frequency distri- 



bution. While this distribution (see Davis et al. 2009 Fig. 14a) is 



biased towards shorter flow lengths as longer jets may become 
invisible in H2 emission once they exit their molecular cloud, 
Davis et aL] (2009) still suggests that the H2 flows in Orion A 



are preferentially short because inclination alone would result in 
a frequency distribution increasing towards longer jet length if 
one assumes that all jets have equal length and are inclined ran- 
domly. |Iomniir&^oebnch] (J20T2a]b| detected 134 MHOs in 
Serpens and Aquila region using the UKIRT telescope and asso- 
ciated them with 131 molecular hydrogen outflows. They found 
that the projected jet length drops exponentially in number for 
longer jets, and does not behave as a power law. They fitted 
the distribution of jet lengths with the exponential function of 
lQ-0J5xiengih{pc) ^^^ ^j^^ fom^j jhat a simple Monte Carlo-type 
model of jets with speeds of 40-130 km s"' and ages between 
4 and 20 x 10^ yr can reproduce the observed length distribu- 
tion. The distribution of jet lengths obtained by us in Ophiuchus 
also shows a decaying distribution for the jets with length > 0.5'. 
However, our sample is too small for a statistical analysis. 

FigurelSlshows the distribution of jet opening angles (middle 
panel) and the distribution of jet position angles (bottom panel). 
Note that the jet position angles show a nearly homogeneous 
distribution (a Kolmogorov-Smirnov (KS) test shows that there 
is a >97% probability that such a distribution is drawn from 
a homogeneously distributed sample), indicative of molecular 
hydrogen outflows in Ophiuchus being oriented randomly. This 
is in agreement with the findings for Perseus ( Davis et al.|2008| l 
and Orion A ( jDavis et al.||2 009) but dilTerent from the results 
from DR21AV75 lDavis et al.,,2007) . [Davis et al.[ ( |2007] l have 
found that the H2 outflows, in particular from massive cores, are 
preferentially orthogonal to the molecular ridge. The alignment 
of position angles of outflows to some preferred direction seems 
to indicate a physical connection between PAs of outflows and 
the large-scale cloud structure ( [Davis et al.|2007) l. 

For the random orientation of outflows in Orion A, [Davis[ 
[et al.[ ( [2009] ) discussed its possible explanation. The orientation 
of outflows could be connected to the position angles of cores 
based on the standard paradigm that clouds collapse along field 
lines to form elongated, clumpy filaments from which chains of 



2009|. To test this paradigm in detail, [Davis et al.| ( [2009| l se 
lected a sub-sample of 22 H2 jets that are associated with 1 . 1 mm 
dust cores. They surprisingly found that even from this finely- 
tuned sample, the distribution of H2 jet position angles with re- 
spect to core PAs is completely random. They suggested that the 
random orientation of outflows may be caused by the relatively 
poor spatial resolution of the millimeter observations which can 
not disentangle the molecular envelopes associated with multi- 
ple sources. 

Of our 21 jets in Ophiuchus, only 4 jets are associated with 
1.1 mm dust cores. We checked these four jets and found that 
there is no obvious correlation between the orientation of out- 
flows and the position angles of dust cores, which agrees with 
the results from [Davis et al.[ ( [2009l l. 

Figure l6] (top panel) shows the relation of jet lengths vs. 
spectral indices of the driving sources, which does not exhibit 
any significant correlation. The spectral index of the longest jet 
has a value of -1.35. For the Class source VLA 1623-243, 
which drives the second longest outflow with a length of 6.2', 
a spectral index from fitting the SED from 2jjm to 24jum can 
not be obtained as this source is not detected in the IRAC ob- 
servations. Of course, we cannot expect a simple linear rela- 
tion between jet lengths and spectral indices of driving sources 
because the relation between them should be complex: for the 
very young sources, we should expect a short jet length as a jet 
should start from the length of zero; then jets probably expand, 
and later on get shorter again as force support fades (Stankd 
2000). However, we didn't see any similar correlation between 
jet lengths and spectral indices of driving sources from Figure|6] 
either. The mean length of jets driven by protostars (a > -0.3, 
including Class O/I and Flat sources) is about 2.1', which is ba- 
sically the same as the mean length of jets driven by Class II/III 
sources of about 2.1'. Similar results have been found in Perseus 
( [Davis et al.|2008) l and Orion A ( [Davis et al.|2009l l. The H2 jets 
in Perseus and Orion A also show no correlation between jet 
lengths and flow source spectral indices. 

Millimeter observations indicate that protostars, which are 
in an earlier stage of evolution, have an increased likelihood to 
drive more collimated CO outflows ( |Lee et aL][2002[ [Arce &[ 
Sargent"2006). For molecular hydrogen outflows, however, no 
such trend has been observed. [Davis et al.[ ( [2008]|2009i) found no 
correlation between jet opening angle and the age (via the proxy 
of spectral index) of the driving source in Perseus and Orion 
A. Our sample of molecular hydrogen outflows in Ophiuchus 
also does not exhibit a correlation between jet opening angle 
and flow source spectral index (see Figure |6] bottom panel). In 
Figurel6] eleven previously known H2 outflows, which have both 
jet length and jet opening angle measurements, are marked with 
black filled circles. Five previously known H2 outflows, which 
have only one H2 feature for each and thus their outflow open- 
ing angles can not be calculated, are labeled with black open 
circles. The three outflows that consist of newly detected MHOs 
are marked with red filled squares. These three outflows also 
each have only one H2 feature. Note that in the top panel there 
are 3 outflows with the same jet length (I'l) and similar spectral 
indices (-1.26, -1.2, and -1.19, see table. [3]l. Note also that there 
are two outflows without spectral indices because their driving 
sources are invisible in the IRAC bands and thus their spectral 
indices can not be obtained by fitting the fluxes from l/im to 
24yum. For the molecular hydrogen outflows in Ophiuchus, the 
mean opening angle of jets driven by protostars (a > -0.3) is 
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about 15.7 degs and the mean opening angle of jets driven by 
Class llfUl sources (a < -0.3) is about 13.6 degs. 

Davis et al. ( 2009|l analyzed the reasons for a lack of cor- 



Jel Length (ofcrhin) 
4 S 



relation between a and jet length or jet opening angle. Firstly, 
shock-excited H2 emission is not a good tracer for outflow pa- 
rameters because of its short cooling time, making H2 emission 
not a sensitive tracer for long jets. In addition, the wings of jet- 
driven bow shocks are often much wider than the underlying 
jet and changes in flow direction due to precession. These facts 
could result in the under-estimation of jet length and the over- 
estimation of jet opening angle. Secondly, the precise relation- 
ship between source spectral index and source age has been not 
established, yet. Our sample in Ophiuchus is somewhat biased 
towards shorter jets compared to Perseus or Orion A because 
our observations just cover part of the p Ophiuchi cloud. On the 
one hand, our survey only covers part of area of the L1688 core, 
while the survey of DavisetaLI (12008) covers the whole western 
portion of Perseus and the survey of Davis et al. ( 2009) covers 
almost the entire Orion A molecular ridge; on the other hand, the 
typical field of view of Sofl is about 0.2x0.2 pc^ at the distance 
of 119 pc, while the typical field of view of WFCAM is about 
6x6 pc^ at the distance of 420pc and about 5x5 pc- at the dis- 
tance of 300 pc. Therefore it is more unlikely for our sample to 
exhibit a correlation between a and jet length. 
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Fig. 4. Distribution of outflows in Ophiuchus. The driving 
sources of outflows identified in this paper are marked with red 
filled pentagrams whose sizes are proportional to the length of 
the outflow. The YSOs identified by 'Evans et al. ( 2009 1 are la- 



beled with open gray circles and the millimeter sources identi- 
fied by |Young et al.| ([2006) are marked with blue pluses. The 
open gray squares show the distribution of known Herbig-Haro 
objects (from the SIMBAD database). The region that we have 
selected for outflow statistical study is outlined with a black solid 
polygon. Note that for the YSOs and millimeter sources only 
those that are located inside this polygon are plotted. 



4. Discussion 

4.1. What kind of YSOs drive H2 outflows? 

Of our 2 1 driving sources of the molecular hydrogen outflows, 7 
(33%) are Class 0/1 sources; 2 (10%) are flat spectrum sources; 
1 1 (52%) are Class 11 sources and 1 (5%) is a Class 111 source 
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Fig. 5. Distribution of jet lengths (in 0.5' bins) (top), jet open- 
ing angles (in 2.5° bins) (middle), and jet position angles (in 20° 
bins) (bottom). The open columns are for all molecular hydro- 
gen outflows in Ophiuchus identified by this paper and the gray 
filled columns are for the previously known molecular hydrogen 
outflows in Ophiuchus. 



(criterion for YSO classification is from 'Gre ene et al.||1994) l. 
FigurelTlshows the distribution of spectral indices for the outflow 
sources and all the YSOs in the "complete" region shown in Fig. 
fflwith a solid polygon. Obviously, outflow sources and all YSOs 
have a similar distribution of spectral indices (a Kolmogorov- 
Smirnov (KS) test shows that there is a >95% probability that 
these two distributions are drawn from the same distributed sam- 
ple), which means that in Ophiuchus the molecular hydrogen 
outflows can be driven by YSOs at any evolutionary stage. This 
result is different from the previous studies. 

In Perseus, Orion A and other distant high mass star forming 
regions such as W75/DR21, the majority of molecular hydro- 
gen outflows are driven by protostars with the positive spectral 
indices ( [Davis et al . 2007; Kumar et al. 2007; Davis et al.|2008| 



2009 1. The mean values of spectral indices of outflow sources for 
Perseus, Orion A and W75/DR21 are 1.4, 0.86 and 1.9, respec- 
tively. The mean value of our flow source sample in Ophiuchus is 
-0.52, which corresponds to Class 11 sources. However, there are 
several uncertainties in determining YSO spectral indices. The 
YSO spectral indices in WR75/DR21 were obtained using only 
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Fig. 6. H2 jet length (top) and jet opening angle (bottom) plot- 
ted against source spectral index, a. The black filled circles rep- 
resent the previously known H2 outflows which have both jet 
length and jet opening angle measurements. The black open cir- 
cles represent the previously known H2 outflows which have 
only jet length measurements. The outflows that consist of newly 
detected MHOs are marked with red filled squares. 



IRAC photometry, which would result in the over-estimation of 
the mean value of a (Kumar et al.||2007[ [Davis et aL]|2009| l. 
Moreover, the source spectral indices used in this work are the 
de-reddened values from [Evans et al.| ( |2009[ ). In fact, if we use 
the apparent spectral indices, the mean value of a of flow sources 
in Ophiuchus will change to -0.16. Certainly, even though con- 
sidering these uncertainties, the value of -0.16 is still signifi- 
cantly lower than the mean values observed in three other star- 
forming regions. Actually, the biggest difference between our 
flow source sample in Ophiuchus and other samples in Perseus, 
Orion A or W75/DR21 is that more than half of the 21 driving 
sources of H2 jets in Ophiuchus are Class II sources. 

Results from the Spitzer c2d project have shown that the 
Ophiuchus and Perseus molecular clouds have nearly the same 
proportion of protostars (Class I and flat-spectrum sources) to 
Class II sources, i.e. 40% in Ophiuchus and 45% in Perseus 
(see Table 5 in [Evans et al.|20()9] l. Therefore, the observed dif- 
ference between Ophiuchus and Perseus regarding the associa- 
tion of molecular hydrogen outflows with protostars and Class 
II sources is not caused by different proportions between proto- 
stars and Class 11 sources in these two regions. The difference 



between the Ophiuchus molecular cloud, where most molecu- 
lar hydrogen outflows are driven by Class II sources, and the 
Perseus, Orion A, and other distant high mass star forming re- 
gions such as W75/DR21, where the majority of molecular hy- 
drogen outflows are observed to be driven by protostars, can be 
explained with the observational selection effects. 

We converted the fluxes of the MHO features in Ophiuchus 
to the H2 1-0 S(l) line luminosity with the convers ion of Lh^ = 
Fh2 X ^ncP-- We adopted 119pc ( [Lombardi et al.| [2b08 ) as the 
distance {d) of Ophiuchus. Note that we did not consider the ef- 
fect of extinction. Figure. [8] shows the relation of H2 1-0 S(l) 
line luminosities vs. spectral indices of their driving sources. 
We found that the median of H2 1-0 S(l) line luminosity of 
the MHO features in Ophiuchus is about 3x10 *■ Lq. Note that 
there are two YSOs ([EDJ2009] 800 and VLA 1623-243, see 
Table. [3]l with no spectral indices because these two sources are 
invisible in IRAC bands and their spectral indices can not be ob- 
tained by fitting the flux from 2//m to 24yum. Here we give the 
value of 2 for the spectral indices of [EDJ2009] 800 and VLA 
1623-243 just in order to plot the H2 line luminosities of the 
MHO features driven by them on the figure. The lines with dif- 
ferent colors and styles show the different detection limits from 
[Davis etal.[ ( |2009ll, [Davis et al.[(|2007] l, [Ioannidis & Froebrich| 
( [2012a|b| l, and [Davis et al.| ( [20()8| l."TTie black sofid fine shows 
the detection limit of the H2 survey towards the Orion A molec- 
ular ridge d Davis et al.| 2009 ). The dot-dashed line shows the 
detection limit of the H2 survey towards Perseus ( [Davis et al.| 
2OO8 1 private communication), while the red dotted line shows 



the detection limit of the H2 survey towards DR21/W75 (Davis] 
|et al.[[20(T7l l. The blue dashed line shows the detection limit of 
the H2 survey towards the Serpens/ Aquila region on the Galactic 
plane ( [loannidis & Froebrich|2012a|b| l. Note that we have con- 
verted the flux detection limits to the luminosity units assuming 
that Orion A is located at the distance of 420 pc and Perseus 
is located at the distance of 300 pc ( |Bally et aL]|2008[ ) as well 
as DR21/W75 and Serpens/Aquila are both located at the dis- 
tance of ~3 kpc. Obviously, most of MHO features in Ophiuchus 
have lower H2 1-0 S(l) line luminosities than the detection lim- 
its of the surveys towards Orion A, Perseus, DR21/W75, and 
Serpens/Aquila. In other words, the p Ophiuchi molecular cloud 
is the nearest star forming region to us and therefore, we have a 
larger possibility to detect the faint molecular hydrogen outflows 
from T Tauri stars. With deeper H2 imaging we can expect to de- 
tect many faint MHOs driven by T Tauri stars in Perseus, Orion 
A and other distant star forming regions such as DR21. 

There are 123 YSOs in the L1688 core (just our "complete" 
region which is shown in Fig.|4]i, including 19 Class O/I sources, 
30 Flat-spectrum sources, 64 Class II sources and 10 Class III 
sources (only the identifications from Spitzer, Evans et al. 2009[ l. 
This classification is based on the de-reddened spectral indices 
of the YSOs in Ophiuchus. In deriving these de-reddened spec- 
tral indices, [Evans et aL] ( [2009[ ) assumed a value of Ay - 9.76 
mag for extinction toward all YSOs in Ophiuchus. So the num- 
bers of YSOs in different evolutionary stages (Class I, Flat, II, 
and III) would change if we adopt another extinction assump- 
tion. 

The value of Ay - 9.76 mag is the mean extinction of the 
whole p Ophiuchi molecular cloud, but our survey only covers 
the partial region of L1688. In fact, the L1688 core has much 
higher extinction than other parts of the p Ophiuchi cloud. The 
extinction map taken from c2d final delivery datc|jshows that the 



'^ http://data.spitzer.caltech.edu/popular/c2d/20071101_enhanced_vl/ 
oph/extinction_maps/OPH_270asec_Av.fits 
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L1688 core region has visual extinction ranging from 10 mag to 
35 mag. For L1688 core, 'Evans et al. (2009i recalculated the 
de-reddened spectral indices for Class I and Flat sources using 
different extinction corrections and they found that about 36% 
flat-spectrum sources in Ophiuchus would de-redden to Class II 
sources if a more realistic extinction correction is applied to the 
estimate of the de-reddened spectral index. Using this percent- 
age for the correction of the number of flat-spectrum sources, we 
estimate that about 11 of 30 flat-spectrum sources are highly red- 
dened Class II sources. After this correction, there are 19 Class 
O/I sources, 19 Flat-spectrum sources and 75 Class II sources in 
the L1688 core region. 

Of course, the discussion above assumes that the same 
fraction of flat-spectrum sources which drives jets as the flat- 
spectrum sources which does not drive jets deredden to Class 
II sources. However, jets may highlight the "true" flat-spectrum 
sources in fact, and only the flat-spectrum sources which do not 
drive jets might deredden to Class II sources. Considering this, 
if we assume that 2 flat-spectrum sources which drive jets are 
"true" flat-spectrum sources and about 10 of 28 flat-spectrum 
sources are highly reddened Class II sources, then there are 
19 Class O/I sources, 20 flat-spectrum sources and 74 Class II 
sources in the L1688 core region. 

Of our 21 outflows, 7 are driven by Class O/I sources, 2 are 
driven by "flat spectrum" sources and 1 1 are driven by Class II 
sources. Now we can deduce that the percentage for protostars 
(Class Qfl+Flat) and T Tauri stars (Class 11) to drive molecular 
hydrogen outflows is ~23% and ~15%, respectively. We can see 
that there is no significant difference between the percentages of 
protostars and T Tauri stars which drive the molecular hydro- 
gen outflows. This might mean that the occurrence of molecular 
hydrogen outflows is not correlated with the age of the driving 
sources. 

Here we should mention the works of |Hatchell et al.| ( [2007| l. 
Hatchell et al. ( 2007| l investigated the CO molecular outflows 
in '^CO (J - 3-2) towards 5 1 submilimetre cores in the Perseus 
molecular cloud using James Clerk Maxwell Telescope (JCMT). 
They calculated the outflow momentum fluxes for all the possi- 
ble outflow sources and they found that outflow power may not 
show a simple decline between the Class to Class I stages. 
These results might indicate that the evolution of outflows is not 
a simple process. 

4.2. Velocity variations along the outflows 

Proper motion studies for HH 34 have found a systematic de- 
crease in proper motions with distance from the driving source 
([Devine e t al. 1997). To explain this velocity decrease, [Cabrit &| 
Raga '{ 2000) modeled HH 34 and found that the variation in ve- 
locity in HH 34 is most likely to be the result of environmental 
drag on the propagation of individual jet knots. 

We also investigated the possible velocity variation along 
outflows in our MHO sample. Out of 21 molecular hydrogen 
outflows, 5 single-lobe outflows have at least two features with 
the reliable PM measurements while 2 bipolar outflows have at 
least two features with the reliable PM measurements in any one 
of two lobes. For each lobe of these 7 outflows, we calculated the 
standard deviation of the velocities. We found that the median 
value of the standard deviations is about 8 km s"' and the mini- 
mum and maximum values are 4 and 47 km s"', respectively. In 
order to describe the variation in velocity along outflows (VV), 
we fitted the distribution of velocities of MHO features with the 
distances from the driving sources using a linear least square 
method. For the bipolar outflows such as the outflow driven by 
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Fig. 7. Histogram of the de-reddened spectral indices (in bins 
of 0.3) of aU YSOs ( [Evans et al.|2009| l located inside the poly- 
gon on Fig.l4](open columns), all H2 flow sources identified by 
this paper (filled gray columns), and H2 flow sources excluding 
those associated with the three newly detected flows (filled black 
columns) . The dashed lines show the criteria of YSO classifica- 
tion from |Greene et al.|fl994[ ). 
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Fig. 8. H2 1-0 S(l) line luminosity plotted against source spectral 
index, a, for the MHO features in Ophiuchus. The black solid 
line, the red dotted line, the blue dashed line, and the green dot- 
dashed line show the detection limit from [Davis et al.| ( |2009|l, 
Davis et al.'(2007jl, jloannidis & Froebrich[ ( |2012a|b| l, and |Davis[ 
et al. (2008, private communication), respectively. 



VLA 1623-243 and the outflow driven by YLW 52, we calcu- 
lated the VVs for each lobe separately. The negative value of 
V V means proper motions of MHO features decreasing with the 
distances from the driving sources. Fig.l9]shows the variation in 
velocity along outflow plotted against the de-reddened spectral 
index of outflow sources. The VVs of each lobe for the bipolar 
outflows are marked with the red circles while the VVs for the 
single-lobe outflows are labeled with the black triangles. The er- 
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ror bars represent the uncertainties of the linear fittings. We also 
exclude a value of VV with very large error bar (the outflow 
driven by [GY92] 93). Note that for the bipolar outflow which 
consists of MHO 2105 and MHO 2106, its driving source is a 
Class source (VLA 1623-243) whose spectral index can not be 
obtained due to the lack of IRAC detection. Here we give a value 
of 2 for the spectral index of VLA 1623-243 with the aim to plot 
the VVs of the outflow driven by VLA 1623-243 in the figure. 

Given the fitting error, two H2 jets show significant decrease 
in proper motions along the outflow. Most of molecular hydro- 
gen outflows in Ophiuchus have roughly constant proper mo- 
tions with increasing distance from the outflow sources. This 
result agrees with the finding by McGroarty et al. ( 2007| l. 
McGroarty et al. ( 2007| l investigated the proper motions of sev- 
eral HH objects from classic T Tauri stars (CTTS) and they 
found that the velocity appears to be constant despite the large 
distances between consecutive HH objects. The maximum dy- 



namic age of the HH objec ts investigated by McGroarty et al. 
( |2007t is about -10^ years. [McGroarty et al!] ( |2007 |l suggested 
two possibilities for this result: one is that the ejection velocity 
at the source was much higher ~10^ years ago when the more 
distant objects were ejected and has decreased over the interven- 
ing years to its current velocity. Meanwhile the velocity of the 
older/more distant objects has decreased over time via interac- 
tions with the parent cloud; the other is that the velocity at the 
source has remained at approximately the same value over the 
last 10^ years and the distant objects have not been slowed down 
much by their interaction with the ambient medium.The maxi- 
mum dynamic age of our molecular hydrogen outflows is about 
lO"* years. Based on our results, we also can not confirm or ex- 
clude any one of these two possibilities. 
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Fig. 9. Variation in velocity along outflow plotted against de- 
reddened spectral index of outflow source. The VVs of each lobe 
for the bipolar outflows are marked with the red circles while the 
VVs for the single-lobe outflows are labeled with the black tri- 
angles. The error bars represent the uncertainties of the linear 
fittings. 



4.3. Future work 

The p Ophiuchi molecular cloud is the nearest star forming re- 
gion. It is an ideal laboratory for the studies of star formation. We 
identified 21 driving sources for 27 MHOs in Ophiuchus, which 
offered a good sample for future high-resolution observations. 
According to jet launching models, the jet launching zones of 
MHOs are within 0.5 AU from the driving sources, which cor- 
respond to around 40 mas at the distance of Ophiuchus. This 
requirement for angular resolution can be achieved, e.g., with 
ground-based interferometric instruments like GRAVITY, which 
is an adaptive optics assisted, near-infrared VLTI instrument for 
precision narrow-angle astrometry and interferometric phase ref- 
erenced imaging of faint objects ( jEisenhauer et al.|2008) l. 

With an accuracy of 10 micro-arcsecond for astrometry and 
4 mas for resolution, GRAVITY/VLTI could study in detail the 
environment within 0.5 AU from young stars. The largest proper 
motion of MHO features which we detected in Ophiuchus is 
about 247 mas/yr , which corresponds to about 28 yuas per hour. 
Therefore for some MHO features with high proper motions 
we can expect to study the jets "in real time" with the resolu- 
tion of GRAVITY/VLTI. The molecular hydrogen outflows dis- 
cussed in the present paper could serve as targets for high preci- 
sion astrometric monitoring in order to study the launching zone 
of MHOs, the collimation mechanism for central jets, and the 
physics of shocks. 



5. Summary 

We performed a deep near-infrared search for molecular hy- 
drogen emission line objects (MHOs) towards the p Ophiuchi 
molecular cloud, covering an area of ~ 0.11 deg^. In total, we 
discovered six new MHOs and detected 32 known MHOs. Using 
previously-published H2 images, we measured the proper mo- 
tions for the H2 emission features in 32 MHOs. The proper mo- 
tions lie in the range of 14 to 247 mas/yr, corresponding to the 
transversal velocities of 8 to 140km/s with the median tangential 
velocity of about 35 km/s. 

Based on outflow morphology and PM analysis, we associ- 
ated 27 MHOs with 21 driving sources. Out of the 21 outflow 
sources, two also drive CO outflows, four (19%) have 850 //m 
and four (19%) have 1.1 mm dust core counterparts. The molec- 
ular hydrogen outflows have the median length of ~0.04pc and 
are oriented randomly. H2 jet lengths (L) and H2 jet opening an- 
gles (6) show no correlations with the spectral indices a of the 
driving sources, which serve as a proxy for the ages of the out- 
flow sources. 

Among the 21 molecular hydrogen outflows identified in 
Ophiuchus, 11 outflows emanate from Class II sources. Taking 
the high extinction in our H2 emission survey region into ac- 
count, we estimate that 23% and 15% of the protostars and 
T Tauri stars, respectively, drive molecular hydrogen outflows. 
The comparable percentages for protostars and T Tauri stars in 
Ophiuchus which drive H2 outflows are in contrast with previous 
results in the Perseus, Orion A, and W75/DR21 clouds, where 
most of the H2 outflows are driven by protostars. The difference 
can be explained with an observational selection effect, as most 
of the MHO features detected by us in Ophiuchus have low H2 
1-0 S(l) line luminosities which are under the detection limit 
of the surveys towards Orion A, Perseus, DR21/W75, and the 
Serpens/Aquila region on the Galactic plane. 

In most cases the proper motions of MHO features in molec- 
ular hydrogen outflows do not vary with distance from the out- 
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flow source, which agrees with the previous studies such as 



McGroartyetal.(2007i. 
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Appendix A: Description of the 1-12 outflows in 
Opiiiuchus 

We detected six new MHOs and 32 known MHOs, with 107 
MHO features in total. We obtained the proper motion measure- 
ments for 86 MHO features. Table IaTI lists the information of 
each MHO feature. 

In this section, we give a brief description of each MHO. 
Here we offer typically four images for each MHO, one H2 im- 
age with the MHO and YSO identifications, one H2 image with 
the proper motion measurements marked, one Ks image for com- 
parison and one continuum-subtracted H2 image to emphasize 
the MHO identifications. 

- MHO 2102 : The region of MHO 2102 is shown in fig- 
ure |AT] This MHO consist s of tw o knots, MHO 2102a and 
MHO 2102b. pomez et al.| ( |2003) firstly identified this ob- 
ject in their H2 narrowband images as [GSWC2003] 19 and 
pChanzady an et al.| ( [2004| associated it with HH 79 and HH 
711 to the north according to the bow-like morphology of 
MHO 2102b. We obtained the proper motion (PM) measure- 
ments of this object. However, due to the low signal-to-noise 
ratio (SNR) of MHO 2102a in the first epoch image, the PM 
of MHO 2102a is uncertain. The PM vectors of MHO 2102 



point back to the YSO [E DJ2009] 800 identified by [Evans 
let al.| ( |2009[ ). |Johnstone et al., ( .2000, ) also detected a sub- 
millimeter source about 7" to the southeast of [EDJ2009] 
800. Based on the measured proper motions, we suggest 
[EDJ2009] 800 as the likely driving source of MHO 2102. 
MHO 2103 : The region of MHO 2103 is s hown in fig- 
ure A. 2 This object was firstly detected by ,Khanzadyan 



eta 



■|( |2004[ ), who identified five H2 emission line knots on 



H2 narrowband images. In our Sofl (2007) images, we de- 
tected three H2 emission knots, including one new feature, 
MHO 2103c. MHO 2103a and b correspond to the previ- 
ous detections of [KGS2004] flO-Ole and f, respectively. 
Note that [KGS2004] flO-01a,b are located outside the cov- 
erage of our H2 images and [KGS2004] flO-01c,g are too 
faint to show up in our H2 images. We rejected [KGS2004] 
flO-Old as H2 emission feature because it is invisible on our 
continuum-subtracted H2 image, and we suggest it might ac- 
tually be two barely resolved stars based on their morpholo- 
gies in the Ks image. We did not obtain PM measurements 
for MHO 2103 because of the lack of previous epoch Sofl 
data. Khanzadyan et aL\ ( 2004 1 proposed YLW 3 1 as the 



possible driving source due to the accurate alignment of the 
MHO 2103 knots together with MHO 2104 knots (see the 
description about MHO 2104) with respect to YLW 31. We 
note that MHO 2103b shows a nice bow shock morphology. 
Based on the alignment and morphology of the MHO 2103 
features, we accepted the suggestion from'Khanzadyan et al. 
(2004) that YLW 31 is the driving source of MHO 2103. 
MHO 2104 : Figure |A3] shows MHO 2104 and YSOs in the 



region. This object was also firstly detected by [Khanzadyan 
|et al.| ( |2004| l, who identified two H2 emission line knots on 
H2 narrowband images. In our Sofl (2007) images, we de- 
tected four H2 emission knots, including two new features, 
MHO 2104b and c. MHO 2104a and d correspond to the pre- 
vious detections of [KGS2004] flO-Olh and i, respectively. 
We did not obtain PM measurements for MHO 2104 because 



of the lack of previous epoch Sofl data. Khanzadyan et al. 
( [2004) proposed YLW 3 1 as the possible driving source due 
to the accurate alignment of the MHO 2104 knots with re- 
spect to YLW 31. Based on the alignment and morphology 



of the MHO 2104 features, we agree with the identification 
by [Khanzadyan et al.[ ( [2004| l that MHO 2104 is driven by 
YLW 31. 

MHO 2105 : Figure|A3]shows the region of MHO 2105. We 
detected 27 H2 emission features in MHO 2105 including six 
newly detected features (located in the red boxes 11, 111 and 
IV in figure [A3] |. We also marked the position of the well- 
known Class source VLA 1623 in Figure [A3] which is the 
driving source of a well-coUimated high- velocity CO molec- 
ular outflow ( [Andre et al.|1990[[Dent et al.|1995MNakamura[ 
[et al.[[201l| l. MHO 2105 is assoc iated with the red-shifted 
lobe of the CO molecular outflow. ' Davis & E isloeff'er( |l995| l 
and D ent et al.,(1995) discussed in detail the association be- 
tween the CO molecular outflow and the H2 emission fea- 



tures in this region. Gomez et al. ( 2003 1 suggested VLA 
1623 as the driving source of most H2 features in MHO 
2105. jLucas et al.[ ( 2008^ also detected part of the features 
in MHO 210 5 using flieir UKIR T H2 images (see B and D 
in Fig. 19 of Lucas et al.[[2008 1. In this paper, we obtained 
the PM measurements for 24 features in MHO 2105 and the 
results of the PM analysis support the conclusion of [Gomez 
et al.[ ( [2003) . Figure A. 4 A. 5 and A.6[ show the zoomed in 

Note th at for MHO 
( [2006) 1 derived the 



images for the red boxes in figure 



2105c,d,g,h,o,v, Caratti o Garatti et al 



PMs by comparing their H2 images (March 2003) with those 
published by Davis & Eisloeffel[ ( [T995 ]l (April 1993). Our 
PM measurements agree with their calculations with the ex- 
ception of M HO 2105c. The position angle of MHO 2105c 
calculated by Caratti o Garatti et al. ( 2006|l is about 243 de- 



grees. Caratti o Garatti et al. (2006) suggested that the driv- 
ing source of MHO 2105c might be another YSO located 
northeast of MHO 2105c rather than VLA 1623. However, 
our PM measurement for MHO 2105c based on the higher 
resolution images gives the position angle of 329 degrees, 
which supports the conclusion that VLA 1623 is the driving 
source of MHO 2105. 

MHO 2106 : Fig ure [XT] shows the region of MHO 2106. 
Dent et al.[ ( [1995| l firstly de tected three b right knots in MHO 
2106 in their H2 images. Gomez et al. (2003) identified 



10 features of MHO 2106 using the deeper observations. 



Khanzadyan et al. (2004j| also investigated this region, in- 



creasing the number of the detected features in MHO 2106 
up to 13. Lucas et al. (2008 1 detected some features based on 
UKIRT H2 images (named E in Fig. 19 of their paper). Our 
observations of this region are the deepest to date. We con- 
firmed 14 H2 emission features in MHO 2106. MHO 2106 
is associated with the blue-shifted lobe of the CO molecular 
outflow, which is driven by VLA 1623 (Andre et aL][1990t 



Davis & Eisloeffel[[T995l [Dent et a l.|1995; Nakamurai et aLJ 
201 l| l. [Gomez et al.| ( [2003[ ) and [Khanzadyan et al. (2004| 



both suggested VLA 1623 as the driving source of MHO 
2106 based on the morphology and the association with the 
CO outflow. We obtained the PM measurements for 10 knots 
in MHO 2106. However, the PM measurements for 7 faint 
features in MHO 2106 are uncertain due to the low SNRs for 
these features. The PMs of the three bright knots support the 
conclusion that MHO 2106 is driven by VLA 1623. 
- MHO 2108 : Figu re [AiS] shows the region of MHO 2108. 
Khanzadyan et al.[ ( [2004) detected this object as a bow shape 
structure in their H2 images. They suggested a YSO north- 
east of MHO 2108 as the driving source. Our deeper obser- 
vations reveal a more complicated structure for MHO 2108. 
It consists of a bow-like structure, a faint knot, and a long tail 
connected to the infrared nebular around WL 18. The opti- 
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cal counterpart of MHO 2108 is H H 673 ( |Phelps & Barsoiiyl 
2004) l. iPiieips & Barsonyl p004| also detected HH 673b 



40" southeast of WL 18. They suggested that HH 673 and 
HH 673b belong to a bipolar outflow driven by WL 18. We 
obtained a PM measurement for MHO 2108, but the PM is 
uncertain because of the difference of morphology between 
the images in the two epochs. Based on the morphology 



of MHO 2108, we accepted the assumption from Phelps & 
Barsonyl (2004|l that WL 18 is the possible driving source of 



Figure A. 8 also shows MHO 2109, which ex- 



MHO 2108. 

MHO 2109 

hibits a bow shape in our H2 images, pomez et ar'(2003 ) and 
[Khanzadyan et al. (,2004j bo th detected this object in their 



images. Gomez et al.|(2003jl suggested that WL 15 (Elias 



2-29, located 2'4 to the east from MHO 2109 and situated 
outside of our image) as the driving source of MHO 2109. 
WL 15 is associated with a bipolar CO outflow ( |Bontemps| 
[eraL][T996l [Sekimoto et"aL]|1997 ) and the red -shifted lobe 



roughly aligns with the direction of MHO 2109. Khanzadyan 
etal.|(|2004!) confirmed the bow shape of MHO 2109 and they 



inferred that MHO 2108 and MHO 2109 arise from a flow 
connecting both knots with a flow direction to the north- 
east. |Nakamura et al. (2011 1 investigated the CO outflows 
in L1688 using CO (J=3-2) and CO (J=l-0) mapping ob- 
servations. They identified a CO outflow driven by the YSO 
LFAM 26. MHO 2109 is associated with the red-shifted lobe 
of this CO outflow. We obtained the PM measurement of 
MHO 2109. We found that the PM vector of MHO 2109 
points to the direction of the northwest, which indicates that 
the driving source should be located in the southeast of MHO 
2109. There are some YSOs in the southeast region. Here we 
suggest YLW 5 as the likely driving source just because that 
YLW 5 is the nearest YSO to MHO 2109. 
MHO 2110 : MHO 2110 consists of three knots, MHO 
2110a-c. Figure |A.9| shows the region of MHO 2110. 



Khanzadyan et al. 



( 2004j) detected this object in their H2 im- 



ages. They suggested EM*SR 24 as the driving source of 
MHO 2110 because MHO 2110, together with MHO 2111 
(see the following description) and several HH objects (HH 
224, HH 709 and HH 418, which are situated outside of our 
image), is aligned along a NW-SE direction passing through 
EM*SR 24 (see. [Khanzadyan et al.|2004l Fig. 3). For MHO 
2110, we only obtained the uncertain PMs due to the low 
SNRs of the emission features. Therefore, we cannot asso- 
ciate MHO 2110 with any YSO based on our PM measure- 
ments. 

MHO 2111 : Figure [AlO] shows the region of MHO 2111. 
MHO 2111 consists of two knots: MHO 2111a is a central 
knot with a diff'use nebula. MHO 21 1 lb is an elongated knot. 
MHO 2111 has been detected by optical, near- and mid- 
infrared observations. Its optical counterpart is HH 224S. 



Zhang & Wang ( 2009 ) identified its mid-infrared counterpart 



using Spitzer/IRAC images. Phelps & Barsony](2004i) iden 



tified HH 224S at the location of MHO 2111. They also de- 
tected several HH objects (HH 224N, H H 224NW1 and HH 
224NW2, see [Phelps & Barsony|2004l Fig. 12) to the north- 
west of HH 224S. According to the alignment of the HH 
224 complex, [Phelps & Barsony| ("2004) suggested [GY92] 
193 (located " 
figure 



A.IO 



7'to the northwest and outside the boundary of 
I as the driving source of HH 224S . [Khanzadyari[ 
etal.[ ( [2004 i detected MHO 21 1 1 at H2 2.12jum as the coun- 
terpart of HH 224S, but they did not detect any near-infrared 
counterparts of HH 224N, HH 224NW1, and HH 224NW2. 
Therefore they suggested that MHO 2111 might be driven 



byEM*SR24, as MHO 2111, MHO 21 10 and EM*SR 24 
are aligned roughly along a line. However, Zhang & Wang 
(2009) suggested WSB49, the YSO closest to MHO 2111, 
as the driving source based on MHO 2 1 11 's wide-cavity lobe 
morphology, which faces away from WSB 49. We obtained 
the PM measurements for MHO 2111, found that the PM 
vectors point back to WSB 49. Therefore, we confirmed the 
assumption fromj Zhang & Wang|(2009 ) that WSB 49 is the 
likely driving source of MHO 2111. 
- MHO 2112 : MHO 2112 is a faint knot (see Figure |ATT). 
Khanzadyan et al.[ ( [2004) l detected MHO 2112 in their H2 



images. They suggested that MHO 2112 is driven by YLW 
16 based on the morphology and distribution of MHO fea- 
tures (see |Khanzadyan et al.|2004| Fig. 5). MHO 21 12 is too 
faint to measure the proper motion. We have no enough evi- 
dence to associate MHO 21 12 with the nearby YSOs. 



MHO 2113 ■ MHO 2 1 1 3 is a bright knot (se e Figure [ATT 1 
Gomez et aL]p003) and [Khanzadyan et al.| ( ,2004) both de 



tected MHO 21 13. [Gomez et al.[ ( [2003[ ) suggested that the 
nearest YSO [GY92] 235 as the driving source. Based on the 
morphology and the location of H2 features (see Khanzadyan 
et al. 2004, Fig. 5), [Khanzadyan et al.[ ( |20D4) l inferred that 
MHO 2113 may be part of the large H2 outflow driven 
by YLW 15. We obtained the PM measurement for MHO 
21 13. The PM vector of MHO 2113 indicates that the driv- 
ing source should be located in the north of MHO 2113. 
There are several YSOs in the north of MHO 2113. There 
is no strong evidence to associate MHO 2113 with the cer- 
tain one of them. Therefore, here we accept the suggestion 
from Gomez et ai:(20UT) that [GY92] 235 is the likely driv- 
ing source of MHO 2113. 

MHO 2114 : MH O 21 14 is a knot (see Figure [ATT) . 



Khanzadyan et al. (2004i identified this object in their H2 



images. They proposed an association between MHO 2114 
and the nearby YSO, [GY92] 235. The PM of MHO 2114 
obtained by us indicates that the driving source should be 
located in the north. There are several YSOs in the north 
of MHO 2114, including [GY92] 235 and [GY92] 232. We 
have associated MHO 2113 with [GY92] 235 and it seems 
impossible that MHO 2113 and MHO 2114 constitute a 
bipolar outflow which is driven by [GY92] 235 because 
MHO 2113 and MHO 2114 are located in flie same side of 
[GY92] 235. Therefore, here we suggest another possibility 
that MHO 2114 is driven by [GY92] 232 because the PM 
vector of MHO 21 14 points back to [GY92] 232. 
- MHO 2115 ■ Figure [Al2[ shows the region of MHO 2115. 
[Gomez et aL] (p003) detected six H2 feafiires in MHO 2115, 
coiTesponding to MHO 21 15a-e,g. 'Kh anzadyan et al.[ ( [2004| l 
identified an additional feature, MHO 2115f, in this MHO. 



Lucas et al. (2008 1 also detected some features using the 



UKIRT H2 images (named G in Fig. 19 of their paper). 
BBRCG 24 is a Class 11 source ( [Evans et all[2009l l and 
is the possible driving source of a CO molecular outflow 
(Sekimoto et al.[[1997| l. Based on the morphology and dis- 
tribution of MHO 2115 features, [Gomez et al] ( [2003| l and 



Khanzadyan et al. (2004 1 both suggested that these H2 emis- 



sion knots arise from multiple, overlapping flows. The PMs 
of the MHO 2115 features support this conclusion. We infer 
that MHO 21 15b,d,e,f,g are driven by BBRCG 24 based on 
the position angles of the PM vectors and the fact that MHO 
21 15b,d,e,g are associated with the blue lobe of the CO out- 
flow, while MHO 2115f is associated with the red lobe of 
the CO outflow (see [Sekimoto et al.|T997[ Fig. 2). The PM 
vectors of MHO 2115a,c indicate that MHO 2115a,c may 
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be driven by another YSO and their driving source should 
be located in the southwest of MHO 2115a,c. We associate 
MHO 2115a,c with [GY92] 93, which is loca ted ~ 9.5'to 
the southwest (outside the boundary of figure |A.12[ ). HH 
314 is located ~47"to the southwest of [GY92] 93. MHO 
21 15a,c and HH 314 may constitute a bipolar outflow driven 
by [GY92] 93. 

MHO 2116 : MHO 2116 exhibits a bow shape in 
Figure A.13| Based on the morphology of MHO 2116, 
Gomez etal.f l 2003 1 and Khanzadyan et al. ( [2004 ') suggested 



YLW 14 as the possible driving source of MHO 2116. The 
PM vector of MHO 2116 points back to YLW 14, supporting 
this suggestion. YLW 14 also drives a CO molecular outflow 
(see lSekimoto et al.|1997| Fig. 2) and MHO 21 16 is roughly 
associated with the blue-shifted lobe of this CO outflow. 
MHO 2117 : MHO 2117 show s a bow-like structure in 
Figure [Al3| Gome z et al^ ( [2003) suggested that MHO 2117 
and MHO 2116 share the same powering star and thus be- 
long to the same flow. Khanzadyan et al. (2004i also de- 
tected MHO 2117. They suggested that the driving source 
should lie to the southeast of MHO 2117 based on the bow 
shock morphology of MH O 21 17. [Gomez et al. (2003) and 
Khanzadyan et al. (2004]) also proposed WLY 2-48 (outside 



the boundary of figure A.13 1 as a possible driving source as 



this YSO is associated with a bipolar CO outflow ( jBontemps 
et al.||199"6l. However, Nakamura et al. (2011 1 investigated 



the CO outflows in this region. They associated this CO out- 
flow with Elias 2-32, which is located to the northeast of 
MHO 2117, rather than with WLY 2-48. Our PM measure- 
ment for MHO 21 17 is uncertain as the morphology of MHO 
2117 changed somewhat from 2000 to 2007. Moreover, the 
orientation of the PM vector of MHO 21 17 is not in agree- 
ment with its bow-like morphology. Thus we can't asso- 
ciate MHO 2117 with any ambient YSOs. The PM vector 
of MHO 21 17 indicates that the possible driven source could 
be located in its southwest. Therefore, MHO 2116 and MHO 



2117 may have some connections as suggested by Gomez 
|eraL] ( |2053] l. 
- MHO 2118 : MHO 2118 consists of two faint knots MHO 



2118a and MHO 2118b (Figure A.14i. Khanzadyan et al 



(|2004j| suggested YLW 16, the closest YSO, as die driv- 
ing source. We obtained the proper motion of MHO 21 18b, 
but it is uncertain because MHO 21 18b has the low SNR in 
the first epoch image. Using the Owens Valley Millimeter 
Interferometer Terebey et al.| (|1989[) detected a CO outflow 



driven by YLW 16. MHO 2118 is associated with the red- 
shifted lobe of this CO outflow. Therefore, we also suggest 
YLW 16 as the driving source of MHO 2118. 
MHO 2119 : MHO 2119 shows a knot with a faint tail 



in Figure |A.14 Khanzadyan et al. ( 2004 1 detected this ob- 
ject and they suggested that MHO 2119 is related to MHO 
2114. The PM vector of MHO 2119 points back to YLW 16. 



north of YLW 



Terebey et al. ( 1989) detected an extended, blue shifted lobe 



16 in the direction of MHO 21 19 (see Terebey 
etal.|1989| Fig. 2b). It seems that MHO 21 19 and MHO 21 18 



constitute a bipolar outflow driven by YLW 16. Therefore, 
we suggest YLW 16 as the driving source of MHO 2119. 
MHO 2120 : MHO 2120 is faint in Figure |Al4| Because 



there are many MHOs and YSOs in this region, Gomez et al 
( |2003| l suggested several candidates for the driving source 
of MHO 2120, including YLW 15 and YLW 16. Based on 
the morphology and distribution of MHOs in this region 
[Khanzadyan et al.| (2004) inferred that MHO 2120 is part 
of a larger outflow driven by YLW 15. We obtained an un- 



certain proper motion measurement for MHO 2120 because 
of its low SNR in both epoch images. Nevertheless, the ori- 
entation of the PM vector indicates that the driving source 
of MHO 2120 should be locate to the south of MHO 2120. 
YLW 15 might have physical relations with MHO 2120. 



MHO 2121 : MHO 2121 is an elon gated knot in Figure|AT4| 
HH 674 is its optical counterpart (Gomez et al. 1998; Phelps] 



& Barsony|2004). A s there are many YSOs around HH 674, 
Phelps & Barsony] (12004 ) did not offer any definitive iden- 



tification of the driving source. |G6mez et al. (20031 and 
Khanzadyan et al.] ([2004) both detected this object in their 



H2 images. The latter suggested MHO 2121 being part of a 
bipolar H2 outflow driven by YLW 16 (see _Khanzadyan et aT] 
2004 Fig. 5). The PM vector of MHO 2121 indicates that 
the driving source should lie to the southwest of MHO 2121. 
There are many YSOs to the southwest that are located out- 
side our images. Thus we cannot associate MHO 2121 with 
any ambient YSO. 
MHO 2122 : MHO 2122 is a central knot wifli a faint dif- 



fuse nebula in Figure A. 14 Based on the morphology and 



distribution of the MHOs in this region, Khanzadyan et al. 



(2004) suggested MH O 2122 to be part of a bipolar outflow 
driven by YLW 15 (see [Khanzadyan et al.|2004, Fig. 5). The 
PM vector of MHO 2122 indicates the driving source to its 
north. Therefore, we suggest YLW 47, a Class II source iden- 
tified by Evans et al. (2009 1, as the driving source of MHO 
2122. 



MHO 2123 : Figure ]AT4] shows the region of MHO 2123, 
which consists of three features. MHO 2123a is an elon- 
gated kno t, MHO 2123b is a faint nebula, and MHO 2123c 
is a knot. [Gomez et al.] ( ]2003) detected MHO 2 123a in their 



H2 images and named it as [GSWC2003] 7c. Khanzadyan 
]erari ([^04') detected four features in MHO 2123, includ 



ing [KGS2004] f04-05h, which is too fai nt to confirm in our 
H2 images. Based on our reprocessing of Khanzadyan et alj 
( ]2004| ) data, which resulted in a non-detection of [KGS2004] 
f04-05h, we suspect that it might have been misclassified 
in the original analysis. Khanzadyan et al. ( 20041 suggested 
MHO 2 123 as part of a bipolar H2 outflow driven by YLW 16 
(see Khanzadyan et al.|2004[ Fig. 5). We obtained PM mea- 
surements for MHO 2123a,c. MHO 2123b have no enough 
SNR to measure the PM. The PM of MHO 2123a is uncertain 
due to its low SNR in the first epoch images. Our PM mea- 
surement for MHO 2123c indicates that its driving source 
should lie to the northwest. We associate MHO 2123c with 
YLW 47, the closest YSO in this direction. 
MHO 2124 : MHO 2124 consists of two bright nebular fea- 
tures with MHO 2 124b being more di ffuse than MHO 2124a 
(see Figure A. 14 ).]G6mez et al. ( 2003|) and Khanzadyan et al. 



(2004) both detected MHO 21247Khanzad yan et al.] ( ]2004 i 
proposed YLW 16 as the driving source. The PM vectors of 
MHO 2124 indicate that its driving source should be located 
to its southwest. MHO 2124 and MHO 2121 might be part 
of the same outflow. There are some YSOs to the west, lying 
outside our images. We cannot offer a definite identification 
of the driving source of MHO 2124. 

MHO 2125 : MHO 2125 is an elongated knot (Figure [Al5 i. 
Based on the morphology of MHO 2125, which hints at a 



bow shock pointing towards SB or SW, ^Khanzadyan et al. 
( ]2004] ) suggested that MHO 2125 might arise from a NW- 



SE flow driven by YLW 15 or YLW 16. We obtained the PM 
measurement of MHO 2125, but it is uncertain because of 
the low SNR in the first epoch images. This uncertain PM 
vector points back to the du^ection of WLY 2-55, an X-ray 
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source identified by |Pillitteri et al 



YSO, YLW 45, located 



(outside the boundary of figure A. 1 5 



( |2010) l. There is also a 

1.4' to the south of MH O 2125 

([T996ji 



I. Bontemps et al. 



detected a CO molecular outflow from YLW 45, and MHO 
2125 is roughly associated with the red-shifted lobe of this 
CO outflow. 

MHO 2126 : MHO 2126 is a bipolar outflow consisting 
of two knots MHO 2126a and b. Figure A. 16 shows that 
MHO 2126 is close to [GY92] 344, a "Flat spect rum" source 
( Evans et al.|2009 1 with a bipolar nebular. Khanzadyan et al. 
(2004[l suggested this edge-on disk system with a bipolar 



outflow might be the driving source for MHO 2126. Our PM 
measurements of MHO 2126 support this conclusion. 
MHO 2127 : MHO 2127 is a bright knot in Figure [AJ3 
Khanzadyan et al.| ( p004l l found that MHO 2127 is elon 



gated with the semi-major axis pointing towards YLW 16 
(see Fig. A.14[ l. We did not obtain a PM measurement for 
MHO 2127 although it has a point-like structure and enough 
SNR in both epoch images, which means that the proper mo- 
tion of MHO 2127 is very small. We note that MHO 2127 is 
also located to the northeast of YLW 45 (~4' and outside 



the boundary of figure A. 15 1 and that the red-shifted lobe of 
the CO outflow driven by YLW 45 ( [Bontemps et al.||1996| l 
extends roughly in the direction of MHO 2127. 
MHO 2128: MHO 2128 is an extended bow-like stru cture in 
Figure A. 17 Figure A.18|shows the zoom-in image. Grosso 



etaL| ( |2001| ) detected MHO 2128 in H2 2.12 //m. They sug 



gested that MHO 2128 belongs to two different jets ema- 

|Al4ll 



nating from two Class I protostars YLW 15 (see Fig. 



and YLW 52. Khanzad yan et al.| ( |2004| ) also detected MHO 
2128, together with several H2 features to the east of YLW 
52 (see Khanzadyan et al.|2004 Fig. 8). They suggested all 
the features to be driven by YLW 52. |Lucasetal. ( 2008| also 
detected MHO 2128 in the UKIRT H2 images (named I in 
Fig. 19 of their paper). They suggested that MHO 2128 be- 
longs to a large outflow connecting MHO 2128, MHO 2106, 
MHO 2105, and some other features outside the coverage 
of our SofI observations. We identified six knots in MHO 
2128 and obtained their proper motions. MH0 2128a-f show 
the radial PMs (Figure A. 18 1, indicating that they arise from 



the same driving source. The PM vectors indicate an driv- 
ing source to its east. The proper motion results agree with 
the suggestion from Khanzadyan et al. ( 2004| l and contradict 
the suggestions of MHO 2128 being part of a large outflow 
driven by a source located to its west 
MHO 2129 



MHO 2129 is a k not in Figure \KA1\ wi th a 

em- 



zoom-in shown in Figure A. 19 Note that Figure A. 19 



phasizes the identification of MHO 2129 with the inverted- 
gray-scale images. MHO 2129 corresponds to the identifi- 
cation of [KGS2004] f08-01b by Kha nzadyan et al. (2004|l. 
They also identified other two features, [KGS2004] f08- 
01c and d, in the immediate vicinity of YLW 52. However, 
[KGS2004] f08-01c and d are invisible in our continuum- 
subtracted H2 images. We suspect that these two features are 
the mis-identified noise points. We obtained the PM mea- 
surement for MHO 2129, but it is uncertain because of the 
low SNR in both epoch images. Given the location of MHO 
2129, [Khanzadyan et al.| ([2004| proposed YLW 52 as the 



driving source. 
MHO 2130 : 

Figure 



MHO 2130 consists of two 



A. 17 with a zoom-in shown in Figure A. 20 



features in 

Note 



that Figure A.20[emphasizes the identification of MHO 2130 



with the inverted-gray-scale images. MHO 2130a is an elon 
gated knot and MHO 2130b is a faint nebula. iChanzadyan 



[et al.[ ( [2004l l suggested that MHO 2130 arises from the walls 
of a cavity, possibly excited by oblique shocks in an S-shape 
flow from YLW 52. MHO 2130b is too faint for a PM mea- 
surement, but the SNR of MHO 2130a is sufficiently high. 
We obtained the PM measurement for MHO 2130a and the 
PM vector indicates that the driving source should be located 
in the southwest of MHO 2130a. Here we accept the sugges- 
tion from Khanzadyan et al.[ ( [2004| that MHO 2130 is driven 
by YLW 52. 



MHO 2131 : MHO 2 1 3 1 is an elongated knot in Figure A.17 



MHO 2131 corresponds to the i dentifi cation of [KGS2004] 
f08-01g by Khanzadyan et al. (2004). They also detected 
a feature near MHO 2131, named [KGS2004] f08-01h, but 
[KGS2004] f08-01h is too faint in our continuum-subtracted 
H2 images to be confirmed. 'Khanzadyan et al. ( 2004| sug- 



gested YLW 52 as the driving source of MHO 2131 and our 
PM measurement for MHO 2131 support this conclusion. It 
seems that MHO 2128-2131 constitute an S-shape H2 out- 
flow from YLW 52. 

MHO 2132 : MHO 2132 con sists of three features in 
Figure A. 21 [Stanke et al.[([2006[) detected the Class source 
MMS 126, located right between MH O 2132b and c MMS 
126 drives a CO molecular outflow ( [Stanke et al.| [2006) 1 
and MHO 2132 is well associated with this CO outflow. 



Therefore, [Khanzadyan et al.[ ( [2004[ ) suggested MMS 126 
as the driving source of MHO 2132. Our PM measurements 
for MHO 2132 support this conclusion. Note that the PM of 
MHO 2132c is uncertain because of its low SNR in the first 
epoch images. Actually, MHO 2132 has a mid-infrared coun- 
terpart in the Spitzer/IRAC images (Zhang & Wang"2009]l 
and MMS 126 is also visible in the IRAC and MIPS images. 
The outflow shows more extended and more obvious struc- 
tures in the IRAC images, including a bipolar nebular around 
MMS 126 and the emission to th e north of MHO 2132c and 
to the south of MHO 2132a (see [Zhang & Wang|2009| Fig. 
5 and 6). 

MHO 2137 : Figure [A3[ shows the region of MHO 2137 
with a zoom-in shown in Figure A. 22 Gomez et al. (2003 1 
detected this object, but they did not identify a driving source 
for this outflow. The PM vectors determined by us indicate 
a bipolar outflow from the Class I source [GY92] 30 (JEvansl 
et al.|2009]l, which shows up as an infrared nebula in the Ks 



image. 

MHO 2149 : MHO 2149 is a knot shown in Figure|A3] It is 
newly discovered in our deep H2 images. We note that YLW 
31 is located to the south of MHO 2149, but YLW 31 drives 
a NE-SW H2 outflow (see the description of MHO 2104) 
and it has been confirmed as a single star with a resolution 
of 0.25" ( [Costa et al. 2000; Barsony et al.|2003| l. Therefore, 



MHO 2149 is unlikely to be driven by YLW 31. [Phelps & 
|Bars ony[ ( [2004[ ) detected the candidate HH object Ola within 
about 10" of MHO 2149. They suggested a relation between 
01 and HH 419, which is located ~3' to the southwest of 
MHO 2149. Therefore, MHO 2149 and HH 419 may arise 
from the same outflow. 

MHO 2150 : Figure [A3[ shows the re gion of MHO 2150, 
with a zoom-in shown in Figure A.23[ MHO 2150 consists 
of two features MHO 2150a, a bright knot with a diffuse tail, 
and MHO 2 150b, a dumbbefl-shape knot. Caratti o Garatti[ 
[et al.[ ( |2006) l firstly detected this object, naming the features 
as A and B, but this identification is not included in the 
SIMBAD database. We obtained the PMs of MHO 2150 and 
the PM vectors indicate that MHO 2150a and b constitute a 
bipolar outflow driven by the central source [EDJ2009] 807. 
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- MHO 2151 : Figure |AJJ shows the region of MHO 2151. 
MHO 2151 is a new feature with no previous identification. 
We obtained the PM of MHO 2151 and the PM vector indi- 
cates that the driving source of MHO 2151 should be located 
to the east. There ai' e two YSOs in the east of MHO 2151, 
GSS 32 and LFAM 3. |Kamazaki e"tai:] ( |2003l ) detected a NW- 
SE CO molecular outflow named as AS outflow (see Fig. 5 
in |Kamazaki et alTi|2003 ) which is centered around LFAM 
3. They suggested LFAM 3 as the driven source of the AS 
outflow because LFAM 3 is associated with the 6 cm radio 
continuum emission, which is considered to be an indicator 
of the active ionized jet, and with the 1.3 mm dust contin- 
uum source. However, we found that MHO 2151 is not as- 
sociated with blue-shifted or red-shifted lobe of AS outflow, 
so it seems impossible that MHO 2151 is driven by LFAM 
3. Therefore, here we associate MHO 2151 with the nearby 
YSO GSS 32 based on the alignment of the PM vector and 
GSS 32. We also can see some faint extended structures be- 
tween MHO 2151 and MHO 2102 (driven by [EDJ2009] 
800; see the description about MHO 2102), which may in- 
dicate some interactions between these two H2 outflows. 



MHO 2152 : Figure |A^ shows the region of MHO 2152. 
MHO 2152 is a newly discovered feature with no previous 
identification. We obtained the PM of MHO 2152 and the 
PM vector indicates that the driving source of MHO 2152 
should lie to the southwest. There are several YSOs in the 
southwest of MHO 2152. We have no sufficient evidence to 
associate MHO 2152 with one of them. 



MHO 2153 : Figure |A^ shows the region of MHO 2153. 
MHO 2153 is a newly discovered feature with no previous 
identification. We obtained the PM of MHO 2153 and the 
PM vector indicates that the driving source of MHO 2153 
should lie to the southwest. There are several YSOs in the 
southwest of MHO 2153, including [GY92] 236 and [GY92] 
239. Here we suggest [GY92] 239 as the likely driven source 
of MHO 2153 just because [GY92] 239 is the closest YSO 
to MHO 2153. 



MHO 2154 : Figure A.21 shows the region of MHO 2154, 



which is a faint nebula in our H2 image. Zhang & Wang 



(2009 ) identified a mid-infrared counterpart in the IRAC im- 



ages, named EGO 33 (see |Zhang & Wang,2009t Fig. 5 and 



6). They suggested YLW 58 as the driving source of EGO 
33 based on that MMS 126 has a close association with the 
known outflow detected in the near-infrared and in CO (3-2) 
emission ( |Stanke et al. 2006). They argued that EGO 33 is 
unlikely to be associated with MMS 126 as EGO 33 does not 
coincide with the axis of the known CO outflow. MHO 2154 
is invisible in the first epoch images, thus we could not ob- 
tain its proper motion. Here we accept the suggestion from 



Zhang & WanglpOOg] ) that YLW 58 is flie driving source of 

MHO 2154. 

MHO 2155 : MHO 2155 shows a elongated knot in 



Figure A. 13 MHO 2155 is a newly discovered feature with 
no previous identification. We did not get the PM measure- 
ment for MHO 2155 because it is invisible in the first epoch 
images. It seems that there are some faint structures between 
MHO 2155 and MHO 2117 in the H2 narrowband image, 
which indicates that there may be some physical connections 
between MHO 2155 and MHO 21 17. 
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Fig. A.4. Zoom-in view of box I in Fig. A. 3 Others are the same as Fig. A.l 
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Fig. A.8. Same as Fig. A. 1 but for the region of MHO 2108 and MHO 2109. 
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Fig. A.9. Same as Fig. A. 1 but for the region of MHO 2110. 
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Fig. A.IO. Same as Fig. A. 1 but for the region of MHO 2111. 
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Fig. A.ll. Same as Fig. A.l but for the region of MHO 2112-2114. 
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Fig. A.12. Same as Fig. A.l but for the region of MHO 2115. 
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Fig. A.13. Same as Fig. A.l but for the region of MHO 2116 and MHO 2117. 
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Fig. A.14. Same as Fig. A.l but for the region of MHO 2118-2124. 
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Fig, A.15. Same as Fig. A.l but for the region of MHO 2125 and MHO 2127. 



M. Zhang et al.: Proper motions of molecular hydrogen outflows in the p Oph cloud. Online Material p 22 



44:ao,c'- 




* -1 




» 


- 

* 


Wvm- 




MH02 1 26b 


4S'«),C' 


^ 


% [GY92]344 j 
♦a 


?+''*S-30.(J ■ 


• 


1 

[GY921352 



iS*l/"sq,p' 4B,o' 



Sitjht Ascension 



44'MO' 


► 






P 


■♦■ 


J*'X|.0" 














- 




« 






-?*'«■» Q" 


r*. .- ■ , 


» 


• 




; 1,- 



1 6"2?"'LK' ■; ' 





f?i-ght Asc^nsil?" 




w'oo.o'- 


^^^^^B 


■ 


Mr'lStS/r 


^^^^vl 


H 


ic 


^^^p^MH02 1 26b 


^^^1 


5 


^^^^^^K> 


.^^^^1 


o 


^^^^^^^Btt 


9 


o 


^^H 


1 


-i*"4s'So.O" 


1 


?^H 



"■ao.o" 


Iff 


• 

5Gkm/5 ^ 




» 


« 


4i^jaQ' 


■ 


^Ai 






m[GY92]344 


4 5'oaa" 




vj?;' 


jiUfjoa'' 


m 


# - 


» 


[GY92]352 


n 


^ *. 


.. 1 



i.s'^van.E? iBu' 



Right Ascflraion 



Fig. A.16. Same as Fig. A. 1 but for the region of MHO 2126. 
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Fig. A.17. Same as Fig. A.l but for the region of MHO 2128-2131. 
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Fig. A.21. Same as Fig. A.l but for the region of MHO 2132 and MHO 2154. 
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Fig. A.22. Region of MHO 2137 and zoom-in view of box V in Fig. A. 3 The top panel shows the H2 image with MHOs, YSOs 



and PMs of MHO features being marked. The bottom-left panel shows the Ks continuum image and the bottom-right panel shows 



the continuum-subtracted H2 image with MHO features being marked. Others are the same as Fig. A.l 
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Fig. A.23. Same as Fig. A.22 but for the region of MHO 2150 and box II in Fig. A.3 
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Fig. A.24. Same as Fig. A.l but for the region of MHO 2152 and MHO 2153. 



